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I. INTRODUCTION 



The soft-collinear effective theory (SCET) provides a systematic approach for separating 
hard, soft, and collinear dynamics in processes with energetic quarks and gluons |2J, y, |4( . 
In SCET the infrared physics is described by effective theory fields with well defined mo- 
mentum scaling, which are used to build operators order by order in the power expansion. 
The hard perturbative corrections are contained in the Wilson coefficients which can be 
computed by matching computations order by order in perturbation theory in a s . The sym- 
metries and power counting in the effective theory simplify the derivation of factorization 
theorems and provide a systematic method of treating power suppressed contributions. The 
construction of the complete set of allowed operators for a process is one of the first steps to- 
wards deriving factorization theorems. The operators are constrained by gauge symmetry in 
the effective theory, as well as by heavy quark effective theory (HQET) and SCET reparam- 
eterization invariance (RPI) jH The operators and Wilson coefficients are typically 
coupled by a convolution integral over the large momenta of gauge invariant products of 
collinear fields. In some cases perturbative matching computations are not necessary, since 
RPI gives relations between Wilson coefficients that are valid to all orders in perturbation 
theory. 

Heavy-to-light currents, J = qTb, are important for describing a broad range of processes 
with SCET, including both inclusive semileptonic and radiative decays like B — > X v lv and 
B - X sl laBBa&eeijieH.exclusive semileptonic and radiative decays such 
as B -»• tt£v and B -> -yJ^~E£ujl E3EHHH0, and exclusive hadronic 
decays like B — > inr JiT 




Here we will consider higher order RPI relations for heavy-to-light currents in a theory 
SCETj with ultrasoft (usoft) and collinear fields. Any momentum can be decomposed as 
pfi = n .pn^ j2 + n-pn^/2 + p^_, where n M ,n^ are two light-cone vectors satisfying n 2 = n 2 = 
and n ■ n = 2. The vector n M appears as a label for the collinear quarks £ n and gluons 
A£, and the quantum fluctuations described by these fields are predominately about this 
direction. The collinear modes have momentum scaling as (n ■ p,n ■ p,p±) ~ Q(A 2 , 1, A). 
The usoft modes q us , A^ s have momenta p^ s ~ QX 2 . We also use HQET usoft fields h v for 
heavy quarks, where w M is a velocity label vector with v 2 = 1 (see for example [27], lisl). The 
mass of the heavy quark is denoted by m, Q is a hard energy scale ~ m, and A 1 is the 
SCET expansion parameter. The auxiliary vectors n, n and v break part of the full Lorentz 
symmetry of QCD, and this symmetry is restored order by order in the power counting by 
reparameterization invariance under changes in these parameters. For processes involving 
heavy-to-light currents it is often convenient to work in the special frame where v± = 0, so 
that v M = n-v n M /2 + n-v n M /2 and n-v n-v = 1. 

In HQET it is convenient to formulate the RPI constraints p| to all orders in 1/m by 
constructing RPI invariant operators and then expanding them to generate a chain of related 
operators. These operators start at some fixed order in 1/m, but once the RPI invariant 
form of this operator is known, all higher terms in the chain are determined. The RPI 
symmetries in SCET are richer and typically the constraints are derived order by order 
in A. In this case, higher order operators in the chain are not fully determined until the 
appropriate order in A is considered. 
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Reparamaterization invariance constraints in SCET were first considered by Chay and 
Kim Q. The complete set of SCET RPI transformations were formulated in Lf. and 
used to prove that the leading order (LO) SCET Lagrangian is not renormalized to all 
orders in perturbation theory. RPI constraints on subleading Lagrangians and tree level 
currents to 0(X 2 ) were derived in Ref. [3| (and verified in [3] for a basis with v± ^ 0). 
At O(X), the extension to a complete set of heavy-to-light currents constrained by RPI 
relations including currents that appear beyond tree-level was made in Ref. At this 

order, all Wilson coefficients are constrained by RPI except for one scalar, four vector, and 
six tensor currents, for which the one-loop matching was done in Ref. 29| and independently 
in Ref. J.Sflf . For the currents that survive for v ± = the O(X), RPI relations were verified 
in Ref. [3jJ. To simplify the computation, they considered constraints restricted to the 
projected v ± = surface (from the RPI-* transformation defined later) since this involves 
writing down fewer operators. At 0(X 2 ), the allowed set of field structures for the heavy- 
to-light currents was determined in Ref. jlij]. Four quark operator currents first appear at 
this order. 1 Recently the type-II RPI invariance was extended to include light quark mass 



effects and provide constraints on certain m q dependent operators For heavy-to-light 
currents at 0(X 2 ) a complete basis of Dirac structures and the full set of RPI relations have 
not yet been constructed. 

In constructing subleading operators we combine objects that are individually collinear 
and usoft gauge invariant. The logic which ensures that all subleading operators can be 
organized in terms of these objects relies on the decoupling of usoft gluons from the leading 
order collinear Lagrangian by a field redefinition involving a Wilson line Y 0. In section HH1 
we show that all results are independent of the choice of boundary condition for this Wilson 
line. Processes described by SCET can depend on the path of Wilson lines, but this path is 
determined independent of the choice of boundary condition. 

Our main objective in this paper is to to derive the complete basis of currents at 0(X 2 ) 
by constructing a basis that is valid at any order in perturbation theory and including all 
RPI relations. Results are derived for use in the v ± = frame (and we take m q = in all 
currents). Two combinations of {SCET RPI-I, SCET RPI-II, HQET RPI} are formed which 
leave v± = 0, and are called RPI-* and RPI-$ (section IlICJ) . We call these transformations 
on the surface v± = "projected RPI relations" and study their completeness in the full 
space of allowed transformations in section Hi Dl For the 0(X 2 ) heavy-to-light currents, we 
show that transformations on the projected surface give the complete set of relations for 
currents defined on this surface (see section IHI D|) . By eliminating the field operators we 
show that it is convenient to consider the RPI relations as constraint equations of the form 

J>K)rf = X)Q(^)if , (l) 

i,k j,£ 

where Bi and Tf are Wilson coefficients and Dirac structures for operators that appear at 
some fixed order in A, and Cj and Tf are terms that appeared in operators from lower orders. 
By deriving these constraint equations in section UlI BI prior to searching for their solutions, 



1 In the most common decomposition the Wilson coefficients of the four quark operators start at Q(a 2 s ), so 
these operators are not needed if the basis is restricted to LO in a s (rn&), such as in Ref. 0. 
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it becomes easier to simultaneously consider the restrictions imposed by the five different 
types of RPI invariance from both SCET and HQET, since each gives a separate constraint. 
A simple counting procedure is given for determining all possible Dirac structures prior to 
imposing the RPI conditions. The solution of the constraint equations in section IIII CI give 
relations between the Bi and Cj coefficients and determine the allowed Dirac structures Yf 
in terms of . 



II. INGREDIENTS FROM SCET 



A. Degrees of freedom, power counting, gauge invariance, and Wilson lines 



We briefly review some basic definitions from SCET that we will need for our computa- 
tions. The fields include collinear gluons A%, ultrasoft gluons A£ s , collinear quarks £ n , and 
heavy quarks h v . An important attribute of our collinear fields is that they carry both a 
large label momentum p and a coordinate x, such as £ np (%). The label momenta are picked 
out by momentum operators, V£, n , P = n-p^ p and V±£ n ,p = P±£,n, P (see Ref. p|), while 
derivatives id^ act on x and scale like ultrasoft momenta. We define collinear covariant 
derivatives 

in- D C = V + gn- A n , iD^ = V[+ gA^ , in ■ D c = in ■ d + gn ■ A n , (2) 
and an usoft covariant derivative 

iD*. = id" + gA*.. (3) 
To construct gauge invariant structures, it is useful to define the collinear Wilson line 



W 



ex p 



perms 



9_ 
V 



n ■ A r 



(4) 



and an ultrasoft Wilson line 



Y{x^) = Pexp (ig / ds n-A us (x M + sn 



(5) 



where it is convenient to choose the reference point sq to be so = — oo with P = P path 
ordering for both quarks and antiquarks. We comment omthe Sq independence of results in 
the next section. Making the collinear field redefinitions 



A n , q {x) ^ Y{x)A n , q {x)Y\x) 



(6) 



removes all couplings to usoft gluons from the leading order collinear Lagrangian and induces 
factors of Y in operators (giving a simple statement of usoft-collinear decoupling). 

We will use the following structures, which are both collinear and usoft gauge invariant: 



Xn = W^ n , H v = Y^h v , V c = W^D C W , V us = Y^D US Y 



(7) 
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collinear quark 


soft quarks 


label operators 


covariant derivatives 


Operator 




hi) Qus 


v v*t 


in ■ D c 


iDc^ in ■ D c iDu S 


Scaling 


A 


A 3 A 3 


A A 


A 


A A 2 A 2 



TABLE I: Power counting for effective theory operators. 

as well as the V± label momentum operator. The fields in Eq. ((7|) are all post-field redefi- 
nition. It is convenient to be able to switch the collinear derivatives for field strengths, for 
which we use 



in-T* c = in-d + ign ■ B , 
Here the field strength tensors are 



-V 



igB 



'A' 



in-D c = in- d — ign-B . 



V 



ign ■ B 



V 



[in ■ T> c , in ■ T> c 



(8) 



(9) 



where the label operators and derivatives act only on fields inside the outer square brackets. 
To determine which fields appear in the heavy-to-light current at each order of A, we need 
the A-scaling of the operators listed in Table UJ 

For convenience we will also use the shorthand notation 



Xn,u> — 

(ign-B)^ = 



Xn 5{u-n ■ vV ) 
igB± 5{uo — n ■ vV y ) 
ign ■ B 5(cu — n ■ vV 



(10) 



so that Xn,uj corresponds to the gauge invariant combination of fields (£ n W) carrying large 
O(X ) momentum u. An operator built out of several of these components then has multiple 
labels, J(ui,u>2, ■ ■ •), and the Wilson coefficient for the operator will be a function of the 
same u)i momentum labels, C(co>i,co>2, • • ■)■ 

We will use a T subscript or superscript to denote objects transverse to and a _L to 
denote those perpendicular to n M and n^, 



= R^-v^vR, 



n/' 



n' J 



R» n-R n-R. 

2 2 



(11) 



The effective theory fields satisfy the projection relations P n Xn = Xn, XnPn = Xn, and 
P V H V = TC V where the matrices are 



Pn = L j, P n = L j, and Pv = ^{l + 



(12) 



The number of independent Dirac structures in a current is reduced by these relations. For 
Xn^'Hv, we can project the Dirac structure onto a four dimensional basis {1,7 5 ,7"} using 



T = tr [P n TP v ] + 7 5 tr [^P n TP v ] + ^ tr [^P a rP v 



(13) 
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where = indicates that the relation is true between Xn and 7i v . Similarly, between collinear 
quark fields, Xn^Xn, we can project the Dirac structure onto the basis {yi, -^7 5 , ^7"} using 

T = I tr [iPnYP n ] + ^ tr [^iPnTPn] + tr [^P n TP n ] . (14) 
o o o 

Finally we define = fi p n u & vpo '/2 where e^ 2 = +1, and note that the tensor identity, 

g Q M = -gfe>? , (15) 

will be useful. 

B. Comments on boundary conditions for Y(x) 

It is worth making a few comments on the path and So dependence of the ultrasoft Wilson 
lines used in Eq. (JHJ). This field redefinition is universal and should apply equally well for 
any physical process. We define 

Y"(x M ) = Pexpfig f ds n-A us (x%)\ , (16) 

J so 

Y\x^) = P'exp(-ig j ds n-A us (x%)) , 

Js 

where x p = x M + sn M . With respect to the equation of motion, n-D us Y = 0, the point so 
implements a boundary condition at infinity, and P denotes path ordering P or anti-path 
ordering P. If is to be the hermitian conjugate of Y one requires that Sq = Sq and P = P. 
This ensures that Y^Y = 1 and that the field redefinition in Eq. (JHJ) causes the usoft gluons 
to decouple in the collinear Lagrangian. The following definitions will also be useful 

Y + = Pexp (ig j ds n-A us (x^ , Y~_ = Pexp (-4g jdsn-A^x^ , (17) 

Yl = P exp (-ig j ds n- A us (x%) ) , Yj = P exp (ig j^ds n- A us (x%) ) . 

Here (Y±)^ = K£, and the subscript on Y± should be read as (Y^)± rather than (Y±Y . 
A common choice for sq is the one made in Ref. 0|, 

s = s = -00, P = P, P' = P, (18) 



where Y = Y + and Y^ = Y_. In Ref. |22| the choice so = +00 with P = P was made in 
order to correspond with particle production, Y^ — Y+. A third possible choice is [3] 

■ (<r>\ - ■ Jp = P,P=P for p,pi >0 

s = — 00 sign! V) , s = — 00 signlP 1 ) , < ~ _ ~, _ '_. . (19) 

V ; V ; |P = P,P=P for P,pt <0 v ' 
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t& n-k+ie 




-n-k—ie 



t* n-k—ie 



FIG. 1: Eikonal ie prescriptions for incoming/outgoing quarks and antiquarks and the result that 
reproduces this with an ultrasoft Wilson line and sterile quark field. 



Eq. (|Tnj) still satisfies sj, = s but corresponds to a different choice for particles and antipar- 
ticles. 2 Here Y = Y + , = Y_ for particles, while Y — Y_, Y^ — for antiparticles. To 
see this recall that 

and that if the label momentum is positive n-p > we get the field for particles, £+, and if 
the label is negative n-p < we get the field operator for antiparticles, £~ j^. Although it 
is important to make some choice for s 0) if one is careful then in any physical problem the 
dependence on so cancels. Any path dependence exhibited by a final result can be derived 
independently of the choice of Sq that one makes in the field redefinition. 

Since the dependence on so sometimes causes confusion, we explore some of the subtleties 
in this section, in particular, why it is important to remember that factors of Y, Y^ can 
also be induced in the interpolating fields for incoming and outgoing collinear states, and 
why a common choice for so = Sq is sufficient to properly reproduce the ie prescription in 
perturbative computations. In many processes (examples being color allowed B — > Dn and 
B — > Xg'y) the So dependence of the Wilson lines cancels and the following considerations 
are not crucial. In other processes, however, the path for the Wilson line is important for 
the final result, particularly when these Wilson lines do not entirely cancel. An example of 



this is jet event shapes as discussed in Refs. 



32 



33. 34 



. _ . See also the discussion of path 

dependence in eikonal lines in Refs. |35 . 

First consider the perturbative computation of attachments of usoft gluons to incoming 
and outgoing quark and antiquark lines. The results for the eikonal factors for one gluon 
are summarized in Fig. ^ and can be computed directly with the SCET collinear quark 
Lagrangian (or from an appropriate limit of the QCD propagator). These attachments seem 
to force one to make a particular choice for sq and so, see for example the recent detailed 
study in Ref. ji^. In our notation it is straightforward to show that this choice corresponds 
to 



s 



• (<r>\ - j- ■ JP = P=P, forP,pt> 

- oo sign (P), s = +oosignrP T ) , <~ ~, _ _' . (21) 

^ " =P =P, for V,V* < 



To see this take a quark with label n-p > and an antiquark with label n-p' < 0, and note 



Note that in this case so = — oo sign('P), is an operator. 
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that 



y£n, P = P exp (ig I ds n-A us (x»j)Q p = P exp (ig [ ds n-A us (^))^ p = Y+Q p , (22) 

J — oo J — oo 

L, P Y ] = g p P'exp (-i<? n-4*(a£)) = £+ p Pexp (z<?y n-4««)) = f+ p yj , 
= P exp (ig j ds n-A us (x^£- p , = P exp (-ig j^ds n-A us (x»f)£~ p , = Y_£~ p , , 
^n, P '^ f = CyPexp (-ig j dsn-A us (x»)\ = CyPexp (-ig j dsn-A us {x^)) = Q p ,Yl . 

J — oo J — oo 



This is in ag reement with the Y = YL, = Yl, Y = Y + , = Y$ used in |42j for the 
production and annihilation of antiparticles and the annihilation and production of parti- 
cles respectively. The results in Eq. ()22)1 reproduce the natural choice of having incoming 
quarks/ ant iquarks enter from — oo, while outgoing quarks/antiquarks extend out to +00. 

Although the choice in Eq. (|2Tj) agrees with the ze's in Fig. ^ it causes complications in 
the attachments of usoft gluons to internal collinear propagators. With Eq. ()21|) we have 
Sq 7^ so- Now the field redefinition still induces factors of Y+Y_ = 1 and y!y + = 1 in 
production and annhilation terms in the collinear Lagrangian, but it also induces factors of 
Y+Y + = Yqo and YlY- = Y^ in quark-quark and antiquark-antiquark terms in the action, 
where 

Y 00 = Pexp(igf dsn-A us {x»)\ . (23) 
J —00 

When usoft gluons attach to a collinear propagator with endpoints x and y we must end 
up with a finite Wilson line Y(x, y). In the original collinear Lagrangian (prior to the field 
redefinition) this finite Wilson line is generated by the time ordering of fields in the usoft 
gluon interaction vertices. If a field redefinition is made with boundary conditions satifying 
s o — $0 then the vertices bordering a collinear propagator induce Wilson lines whose Sq 
dependence cancels, leaving this same finite Wilson line. For example, with sq = —00, 
F(-oo,x)F(-oo,0) t = Y{0,x). A choice like that in Eq. (JZQ) is more complicated since 
it violates hermiticity: (£ n )^ = £t prior to the field redefinition, but this is no longer true 
for the and fields after the field redefinition. Correspondingly, the term in the action 
determining the free propagator depends on Y^. Thus, in this case there are Y factors in 
both the propagators and vertices which must be taken into account in order for the path 
ordering not to conflict with the result from time ordering, and give the same finite Wilson 
line. 

Let's adopt the choice in Eq. (|T5j) rather than Eq. (}2~T]) and check that the theory with 
the field redefinition in Eq. (JBJ) still correctly reproduces the results in Fig. ^ for this case. 
Here we have Y = Y + , Y^ = YJ for particles and antiparticles. Thus, the correct ze's are 
obviously reproduced for the incoming collinear lines as well as intermediate propagator 
states. On the other hand, the result for an outgoing quark seems to have the wrong factor 
since comes with a Yj rather than a Y+. However, with the standard definition of an 
outgoing state there is actually an extra Y^ induced by the field redefinition on the out- 
state itself. When we take this factor into account we have Y^Yl = Yl as expected. To see 



S 



this, recall that an outgoing collinear quark state ou t(p\ is generated by a suitably weighted 
integral over (0|£+(o;t), in the large time limit T — > oo for = {T,x). When we make the 
field redefinition this field, ££( x t) generates an usoft Wilson line which extends from our 
reference point sq = — oo to the n-x point for our asymptotic outgoing state (which is +00 
for T — > 00), namely a factor of Y^. A similar argument applies for outgoing antiquark 
states, where we get Y + Y^ = Y_. The same considerations must also be made for hadronic 
bound states where they apply to the interpolating quark/ antiquark fields used along with 
the LSZ formula to define the outgoing state. The factors of Y^ are universal, independent 
of which out-state we choose. There are no additional factors for our incoming states since 
our reference point and T = —00 coincide, Y(— 00, —00) = 1. Once the factors are taken 
into account, the choice in Eq. (fTHJ) correctly reproduces the path for outgoing quark and 
antiquark lines. If we had instead made the choice for s i n Eq. (|T9*j) (which also satisfies 
s o = *o) then we would have Y$ factors for incoming antiquark states and outgoing quark 
states, but the final outcome is the same. Thus the complete result is independent of the s 
choice. 

The above discussion covers usoft interactions from the collinear Lagrangian, but it is also 
worth remarking on the interactions induced by the field redefinition in (possibly non-local) 
operators that are not time ordered. We continue to use Eq. (fTHjl. Here again, the identity 
Y*Y = 1 is important in order to prove the cancellation of usoft gluon attachments. It is 
convenient to adopt a convention where one collects the extra factors of Y^ induced from 
outgoing states together with the Y^s from production fields in these operators. In this 
case if we consider J(x) = f° r the production of a collinear quark and antiquark, 

then instead of writing only the Y_ and Y + from the fields we write J — ► ££Y+ fiY-^n = 
in^in which includes the Y's from any out-state this current could produce. Here the usoft 
interactions in the Y and Y^ lines extend from x to 00 and cancel. For the annihilation of a 
quark and antiquark the lines extend from —00 to x and also cancel, namely YlY + = 1. These 
two cancellations are often sufficient to ensure the decoupling of usoft gluons. For example, 
in exclusive processes we must have color singlet combinations to connect to incoming or 
outgoing collinear hadrons and so we can typically pair up £^ and £J fields in the hard 
scattering operator and make the cancellations manifest. 

If we instead consider an inclusive process like DIS then we have a quark scattered to a 
quark (we consider generic Bjorken x < 1 in the Breit frame). In this case including the Y^ 
from one outgoing quark in the final state gives Ct^ + $ ^+£n where the Wilson 

lines do not seem to cancel. Here in order for the cancellation of usoft gluons to take place it 
is important to either a) take into account all factors of Y^ from the outgoing proton state, 
or b) include the Y^ from one outgoing quark state but note that we are only matching 
cut diagrams for this inclusive process. The choice a) or b) depends on whether we want to 
take the imaginary part at the very end, or from the beginning. For b) the effective theory 
computation has the imaginary part of the hard computation, but the imaginary part also 
effects the collinear operator, where we can denote the cut by a vertical line, |. With our 
initial state for the T-matrix taken on the RHS of the cut, the signs are as in Fig. ^ but on 
the LHS we have the complex conjugate of these expressions, and the above computation 
becomes 

(L)\mn) - (Z n Yl)\@Y+Zn) = Liin- (24) 
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Thus, the usoft gluon interactions also cancel in this case. Alternatively, with a) one must 
keep track of all the lines in the full forward scattering calculation including Y"iP factors 
from all initial and/or final state quarks, and then the Y's in the low energy theor y ag ain 
all cancel. Both ways we arrive at the same final result, (ImC)£ n ^£„ (see Refs. |43[ 44 1 for 
a discussion of DIS in SCET). Similar considerations can be applied to B — > X s ^ in the 
endpoint region. The so dependence cancels, and for this process we are left with a finite 
usoft Wilson line, h v (x)Y(x,0)h v (0). 

To summarize, keeping careful track of the boundary condition so dependence in the 
usoft Wilson line Y, a choice satisfying s = Sq appears to be the most natural (even though 
there will be additional Y^ factors from states). Physical results are independent of the 
choice made for the so reference point. They may still depend on the path of Wilson lines 
in the final result, but this is determined by the universal class of processes described by 
the operator rather than the choice of so in the field redefinition. Similar conclusions hold 
for the path dependence in collinear Wilson lines W. We note that with respect to the 
definitions of the gauge invariant structures made in Eq.(JTUj), the remaining allowed global 
color rotations simply correspond to color rotations at the reference point. We will pick 
the same reference point in W and Y factors. For example, the gauge invariant product of 
fields (Y*h v ) carries a color index in the 3 representation, which by convention is acted on 
by global rotations U(sq), via (Y^h v ) — > U(so)(Y^h v ). These color rotations still connect 
invariant products of collinear and usoft fields. 



C. Reparameterization invariance 

The structure of the currents is constrained by reparameterization invariance, which is 
an invariance that appears due to the ambiguity in the decomposition of momenta in terms 
of basis vectors and in terms of large and small components. The total momentum P^ of a 
heavy quark is decomposed as 

= mQV » + k", 

where mq is the quark's mass, v M is its velocity, and is a residual momentum of order 
mgA 2 . Then the simultaneous shifts 

v"^v^ + ^ and P^P-mQ/^, (25) 

where the infinitesimal (3^ ~ A 2 , can have no physical consequences 0. We refer below to this 
reparameterization invariance as HQET-RPI. The transformation of the field 7i v — > 7i v +57{ v 
induces terms at O(X ) and C(A 2 ), 

s^n v = {imp ■ x)h v , 5 {x ' 2) n v = i- n v . (26) 

There are also reparameterization invariances associated with ambiguities in the decom- 
position of the momenta of collinear fields. Here the total momentum P M of a collinear parti- 
cle is decomposed into the sum of a collinear momentum with (n-p, n-p, p±) ~ Q(A 2 , 1, A), 
and an ultrasoft momentum k^, with (n ■ k,n ■ k, k±) ~ Q(A 2 , A 2 , A 2 ): 

pv = p » + ^ (27) 
= —n -(p+k) + —n -k + (p ± + k±). (28) 
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This decomposition has two types of ambiguity. The first comes from splitting P M into large 
(p) and small (k) components. Thus operators must be invariant under a transformation 
that takes 

V — > V + n ■ £ , in ■ d — > in ■ d — n ■ £ , 

Vl^Vl+ei, id£ -> id" - 1 ± , (29) 

where all operators and derivatives act on one or more collinear fields, and £^ is 0(X 2 ). We 
refer to this reparameterization invariance as SCET RPI-a. Examples of an infinitesimal 
transformation on fields and operators are 

5p Xn = (z£-x) X n, 5i^n = £ a ± , S^P = n-£, (30) 

where n-t = 0. Note that (i£ ■ x) terms only effect ultrasoft derivatives acting on the fields 
since the overall current is evaluated at x = 0. 

The second ambiguity in the decomposition of the momentum of the collinear particles 
comes from choosing the light-cone vectors n and n. An infinitesimal change in these vectors 
which preserves the relations n 2 = 0, n 2 = 0, and n-n = 2, can have no physical consequences. 
The most general infinitesimal transformations of n and n that preserve these conditions 
along with the collinear power counting are 0] 

(I) / n " ~~ * n " ^ (II) / n>1 ~ * n " (III) < n " ~ > ^ °^ n " (31) 



^ n M -> n M {n^—^n^ + e^ { n^ -»• (1 - a) 

where {A^,e^,a} ~ {A 1 , A , A } are five infinitesimal parameters, and n ■ e 1 - = n ■ e 1 - = 
n ■ A ± = n ■ A ± = 0. 

If we start in the frame v± = 0, then transformations (I) or (II) or (HQET-RPI) take us 
out of this frame. A certain combined type I and type II transformation, however, leaves 
v± = |3lj |. We refer to this transformation as RPI-*. We can also form a combined 
HQET and type II transformation that leaves v± = which we refer to as RPI-$. These 
transformations are 




*v - 


- + 




^n A - 


n^ - 


'V (n-v) 2 


v„ - 





W + > (32) 



Jfy -> 



where A 1 - ~ A, j3 ~ A 2 , and j3± is the _L-part of (3t- In defining the $-transformation 
we found that it is more convenient to leave v j_ = by making a transformation on v 
simultaneously with n, rather than simultaneously with n. Under the ^-transformation the 
components of a generic four-vector transform as 

n-V -^nT + A 1 -^ 1 , 

A ± ■ V ± 
n ■ V — > n-V - 



(n-v) 2 



n n \ _ / n» n M 
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To the order we are working we need the following terms from an RPI-* transformation: 



SpVl = ~P, 5f\ign-B) = A ± ■ {i g B ± ) , (34) 



sFHig&D = o, 2 e 0) xn = o 



5f)V» = (n-vfn^A 1 ^, 5^ 



Xn Xr, 



tA 1 



4 ' 

SpiigB^) = (n-v) 2 n^A ± -(tgB ± ), SpV = ~(n-v) 2 A ± ■ V± , 



e 2) Xn = -Xn(igB x + f L ^)^-^ + Xn\ r ^ 5 igB ± -A 

V / V' 2{n-v y l(n-v) z V 

where tit is the transverse part of n, 



71^ = 11^ — n-t)/ = — (n-v) — . (35) 

We will also need the transformation 

5i x2) 5(co - n-vV ] ) = — A ± -V\ 6'(u - n-vV ] ) . (36) 

n-v 

For the RPI-$ transformation at the order we are working we need the following terms: 

5f } H v = (im(3 T ■ x)H v , 5f ] n v = ^H V , (37) 

4 A2) S(lo - n-vV ] ) = -n-Pr P f 8'{u - n-vP^) . 

For the last identity it is straightforward to see that the $-transformation on n does not 
enter until one higher order. We chose to define the RPI-$ transformation to be for v and n 
rather than v and n because of the property that terms with n are often pushed to higher 
order, making the relations derived with RPI-$ more orthogonal to those from RPI-*. For 
example, in order to have a simple form for the <5$r's in Eq. (|T0|) below it is important that it 
is n and not n that transforms. Finally, we note that since all Dirac structures are 0(1), all 
RPI transformations of Dirac structures have the same power counting as the transformation 
parameter, in particular, 5+T ~ O^X 1 ) and 5$T ~ C(A 2 ). 

Finally, note that we will consider the RPI tranformations of all currents prior to making 
the field redefinition in Eq. © so that we do not have to transform Y. However, in order not 
to have to switch our notation back and forth we will write all equations with the operators 
obtained after the field redefinition. This implies that results quoted for the transformation 
of objects involving TC V should be though of as being made for h v , with the field redefinition 
which induces 7i v made only afterwards. 



D. Completeness of Projected RPI 

It is natural to ask if for v± = the transformations RPI-$ and RPI-* in Eq. (|32jl are 
sufficient to give the complete set of constraints that arise from the original SCET type-I, 
SCET type-II, and HQET RPI transformations. The set { RPI-$, RPI-*, SCET-II } forms 
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FIG. 2: Transformation of operators on and off the v ± = surface. Here 0^2 exist for v± = 0, 
while 0-3,4,5 vanish on the v ± = surface. 

an equivalent complete grouping related by linear combinations. To address this question, 
consider splitting all possible operators into two sets, a set {Oj} which do not vanish on the 
v± = surface and a set {Oj} which do. An example is pictured in Fig. El 

Constraints are derived by requiring cancellations among the resulting post- 
transformation set of operators. If we consider an operator Oj then under one of the 
projected RPI transformations, RPI-$ or RPI-*, it transforms into the set {Oj, 0&}. On 
the other hand an operator Oj only transforms back into the set {Oj}. This is a special 
feature of the projected transformations and ensures that relations derived on the v ± = 
surface can not be spoiled by operators which appear away from the surface. It appears that 
we can neglect the Oj operators since they vanish when we project on the v± = plane. 
However it is still possible that we will miss an additional relation between operators on the 
surface, so that the surface analysis will not be complete. 

There are two possible sources that could lead to additional relations beyond those derived 
from projected RPI on the surface. First, under the SCET RPI-II transformation e± ~ A is 
allowed, while in the RPI-* and RPI-$ transformations we only have smaller transformations 
of ft of 0(A X ) and 0(X 2 ). Thus we could miss relations from the more restrictive e± ~ 
A allowed by SCET RPI-II. Note that an SCET RPI-II transformation takes us off the 
projected surface. Second if we project onto v± = then constraints are derived only 
by enforcing cancellations within the set {Oj}. It is possible that an operator O4 exists 
that is obtained from the transformation of two operators 0\ and O2 that are not related 
by transformations on the surface. Enforcing the cancellation of O4 then relates Oi and 
O2. This is pictured by the star in Fig. |21 A related alternative is an operator like O5 
pictured with the box which is obtained from transformations of O12 and O3. If O3 is 
otherwise constrained then this would also constrain Oi^. In cases with multiple operators 
appearing and multiple transformations we must of course consider the linear independence 
of combinations of operators. If an Oj contributes and it is not otherwise constrained then 
this is not of concern, since in the end we discard Oj by projecting onto the v± — surface 
anyway. We will call an operator that vanishes for v± — but that generates a relation 
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between operators on the surface a "supplementary projected operator" (SPO). 3 To check 
for the existence of an SPO we might in general need the full set of v ± 7^ operators. At 
O(X) the comparison of the results derived in Ref. in the full space, to those derived in 
Ref. J3l[ on the surface v± — shows that there are no SPO's at this order. 

For the 0(\ 2 ) heavy-to-light operators considered here we show that there also no SPO's 
in section 1111 PI This is done by a careful choice of our Dirac basis which makes it simpler 
to demonstrate that there are no further type-II RPI relations, and by explicit construction 
for other possible SPO's. Thus, the analysis on the v ± = surface is complete for our 
computation. 

III. HEAVY-TO-LIGHT CURRENTS TO 0{\ 2 ) 

To order A 2 , the operators and Wilson coefficients for the heavy-to-light currents can be 
written as 



where J( kx '(uJi) represents the 0(X k ) terms with dependence on convolution parameters 
u>i. Here the subscript x distinguishes distinct field structures at a given order, and j 
sums over distinct Dirac structures. At 0(A) we know that there are at most two relevant 
convolution parameters i — 1,2, while we will see below that at 0(X 2 ) there are at most 
three. We will consider both scalar, vector, and tensor currents (and the simple extension 
to the pseudoscalar and axial vector cases). When necessary we add an (s), (t>), or (t) 
superscript to the Wilson coefficients in order to distinguish these cases, e.g. B$ . 

We begin in section IHI Al by constructing all consistent field structures for the NNLO 
currents. In section IIII Bl we use reparameterization invariance to derive the constraint 
equations for these currents under different types of RPI invariance on the v ± = surface. 
In section 1111 CI we solve the constraint equations to find the allowed Dirac structures and 
obtain relations among the Wilson coefficients. Finally, in section 1111 Dl we show that the 
results from the v± = surface are equivalent to those obtained if all relations in the full 
space were projected onto this plane. 

A. Current field structures at 0{\ 2 ) 

We first construct a basis of currents that is consistent with gauge invariance and power 
counting and eliminate structures that are redundant by the equations of motion and Bianchi 



3 In the case of type-II transformations, operators like O4 and O5 need not be in the {Oj} class. 



J 



J(0) + J(D + J(2) 



(38) 



^2 jduj Cj(u, m, ix)jf\u, y) + y^l [duJi] B xj (u) h m, ^)jf x \u h /1) 
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UJ 



identity. At LO and NLO the currents are 

J®(u) = Xn,uTH v , (39) 

At NNLO we find that a convenient basis for the set of field structures for the bilinear quark 
operators is 

J {2a) ^) = ^Xn,.Tl ) iVl a H v , (40) 
J {2b \uj) = — Xn^in ■ V us Tt b )H v , 

UJ 

J {2c \u) = -lxn,Jvi sa r? c) H v , 

J (2d \^ = ^Xn^V^U L 



m(uj 1 + uj 2 ) \uj 2 oji J u> 
J (29) M = -i- Xn*n{(ign ■ BU + 2(igB ± U ■ ^l^p } T ( 9 )^ , 

J (2h) ^i,2,s) = 7 \ \ Xn,^gBiU(igB^UT^n v , 
m(uj 2 + LO3) y> 

J (2l) (ui,2,3) = , \ Tr{(igBjU(igB^U] fc^TgW, . 

m{uj 2 -\- UJ3) w 

For a basis of four quark operators we take 

j^Vi^2,^ 3 )= Yl [ 

f=u,d,s 

j(»)( Wll w ai w S )= £ K,. 2 T A T (fcx)X ^ 3 ] [x„^T (Jfcw) Wj (41) 

where the matrices T A are generators of SU(3) with an implied sum on A and xL nas a 
collinear quark with flavor /, whereas Xn carries the flavor of quark from the full theory 
current. We impose the RPI type-Ill invariance in Eq. (|31|) on all operators by multiplying 
by an appropriate power of n-v. The basis in Eqs. (I39I40I41J) is valid whether or not we take 
v± = 0. The v± = choice only effects the basis of Dirac structures. 

The 11 operators in Eqs. (|4()|41jl can be compared with the 15 field structures in the 
basis of Ref . 3| • We have no analog of their j[ 2 ^ 3 7 currents which have an explicit 
because with momentum labels the multipole expansion is performed directly in momentum 
space jiij]. Correspondingly, our and currents have no analogs in their basis. 
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There is a correspondence, j( 2a ' M ) <-> jfj, J^ 2 ^) «-> J< 2 > «_> ^\ 415 , and 



4,6' ^ " J 8,9,10> J " ^13,14,15' 

our j( 2fc ) encodes their j} 2 / and j}g currents. 

In arriving at Eq. (}4T)j) we have used Eq. (jHJ) to switch to a basis with Pj_'s, in-d, and 
field strengths rather than colhnear covariant derivatives in order to give simpler constraints 
from RPI. The basis with covariant derivatives is more natural from the point of view of tree 
level matching and the relation between the two is discussed in section IIVI The prefactors 
in j( 2a_2j ) have been chosen with these relationships in mind, in order to make the matching 
coefficients for the operators simple. The combinations in j( 2 f< 2 9) were chosen because they 
have simpler transformations under RPI. 

Structures were also removed from Eq. (}4T)j) using equations of motion and the Bianchi 
identity. In the effective field theory this gives a valid basis at anyloop order. After 
decoupling the usoft gluons the LO Lagrangian for colhnear quarks is [4| 

4 0) = e n |(^^ c + ^i^t^)^ = x4(^^ c + ^i^)xn, (42) 

so the equation of motion for \ n can be written 

in-dxn = -(ign-B)xn ~ iftp^c Xn , (43) 

where using Eq. (jHJ) the last term can be written as a sum of terms with either two Pj_'s, 
two (igB±ys, or one of each. Eq. shows that a a current x. n in- dHl is redundant by 
the colhnear quark equation of motion and need not be included in the list, explaining why 
we only have and j( 2c \ (Note that in-T> us x n = in-dxn-) As noted in 12j|, this makes 



their current redundant. In j( 2a ) we have restricted the ultrasoft derivative acting on 
h v to be purely transverse since the heavy quark equation of motion is v ■ D us h v =0. 

One can also consider using the colhnear gluon equation of motion. After the field 
redefinition in Eq. (JHJ), the lowest order colhnear gluon Lagrangian is the same as in QCD (if, 
df* = l/(2g 2 ) tr{[iD%, iD v c ]} 2 . Varying df^ + Cg with respect to the colhnear gluon field 
A A ^ and contracting with fi i [T A gives 

xr{0) i 

= n»T A — = -n„ [iD cv , [W£, ilfl] + gT A £ £ . (44) 

"■™-c/ji 9 j 

Next we multiply by on the left and W on the right, use the identity (W^T A W) ®T A = 
T A ® (WT A W*), and label by uo 2 to give 

-9 2 T A Y,[x f nT A MU = {\iV w \m-V c ,m\)^ (45) 
/ 

oj 2 

Multiplying by Xn W on the left and TTL V on the right where Y is some Dirac structure gives 
^ Xn^{wn-BUYH V = -g 2 £ [xU-^T^xU [Xn,«J A TH v ] (46) 
+ ^2 Xn^V^igB^^YHv + 22^3 Xn,^ [(igB±) 

U>3 
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This result can be used to eliminate the current (c^ lj2 ) in terms of and j( 2fc ) if 
desired. We have chosen not to remove this operator since doing so would induce a tree 
level matching contribution for j( 2k \ For listing results it was more convenient to leave all 
four quark operators with coefficients that start at one-loop order, 0(a 2 s ). Since Eq. (}4*6j) 
eliminates a current that will not show up in the constraint equations it does not effect the 
discussion of RPI relations. 

The Bianchi identity in QCD is D^G U(7 -\- D u G afl + D a G ^ = 0. It can be used to eliminate 
terms proportional to igB^ ± = [iDj_,iD±] or 

(47) 

in terms of factors of igB\ or igB\ = [in-D, iDj_]. The Bianchi identity gives [n-D, B^j] = 
[D^_, B v ± ] - [D u ± , B£] so using Eq. © we have 

(igBTJ = ^(WBI) - ^(ig&D + ± [(igB£), [VigBD] - i [(igBl), {VigE®] . (48) 

Thus a heavy-to-light current with (igB^jJ can be matched onto a linear combination of 
j(2/,2e) an( ^ j(2h) w j^] 1 antisymmetric indices in T?R. 

B. Constraint equations from reparameterization invariance 

We derive constraint equations for the allowed subleading currents considering the differ- 
ent types of RPI in turn. 



i 



WHgB^JV 



1. RPI-* at 0(A) 

To set the stage we review the constraints at 0(A) from SCET RPI. To ensure that the 
next-to-leading order current is RPI-* invariant, we must have 

5 f)j(0) +S pj(i) = . (49) 
Computing the various terms in this equation gives 4 

e ) ^ (0) H = W(^r + e ) r)^, 

8pjW>M = Q. (50) 

The terms that must cancel all have a common dependence on Xn,ta> A^, an d T~t v which 
can be factored out. The remaining coefficients and Dirac structures give the constraint 



4 Note the remark on our use of notation at the end of section lll Cl that explains why we do not include the 
transformation of Y. 
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equation: 




(51) 



where the index a is _L, j sums over Dirac structures, and S"T^ is defined through 



^ 1)r o-) = ^ A ^ r o) 



(52) 



2. RPI-% at <D{\ 2 ) 

The only terms in the current whose transformation under RPI-$ leaves uncanceled terms 
are and j( 2a \ We must have 



Now, 



^ 2 )j(0) + ^°)j(2a) =a 

5f ] J (0) M = Xn[~n ■ PtV*5'(uj - n ■ vV ] )]T7i v + X r 



(53) 



'5pT + 



(54) 



Suppressing the common fields Xn,u, an d vector (3^ leads to the constraint equation 




where 



5f% } = #ajr w and 7? = r - fv«. 



(55) 



(56) 



3. SCET RPI-a at <D{\ 2 ) 

The terms in the current that transform under SCET RPI-a are J (0) , J (la) , J (2b) , and J {2c \ 
We must have 

5 f ) J(0) + ^) j(la) + S f) j(26) + 5 (A«) j( 2c) = g (5?) 



Now, 



5i Al )j( la )M = lxn,.e0fa)K, 

^° ) J (2c) H = --wTf c) ^7i l ,. 



(58) 
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This leads to the a constraint equation between 0(X 2 ) and O(X ) 



U) 



(59) 



and a constraint equation between O (A 2 ) and 0(A) 



(aj) 



(60) 



^. SC£T RPI-* at 0{\ 2 ) 
Under RPI-* we must have 

5 (A 2 )j(0) +(5 (A 1 ) J (l) +5 (AO) J (2) =a 

Many of the currents transform under this form of RPI: 



(61) 



(n-v y 



<^±^± _t t igB±-A ± 
■Aj. -V\[n-v) 5'(uj-n-vV ] )^TH v 



+ t[^ ± ^f-\S(u;-n-vV 



2V 





— Xn.ui \ T^a ^ 






^(A°)j(2a,26,2c) 


= 0, 




(n-v) _ / 

Xn,iv I 


e°)j( 2£ )K 2 ) 


1 



m(u;i + uj 2 ) ' 


$(\°) j(2f,2g,2h,2i,2j,2k) _ q 



(62) 

The terms in Eq. (|6^jl can be grouped into two unique field structures, [x_n,w A„ "P^ • • • Hu] and 
[x nj£J A^M • ■ - TCv), which must cancel independently. This gives two constraint equations. 
The terms proportional to Aj^V^ give 



(63) 



E^H( T ?I) + T (l)) = E { - C ^hhlT U) - 2^.(u;)< ? f r a) _ 
i j 



(aj) ; • 



From Eq. (|51jl we know that the index a on ©7^-) must be _L so the last term vanishes. 



Inserting Eq. (JHT|) also simplifies the nonvanishing terms. Finally we know that is 
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symmetric in a and (3. With these simplifications we have the constraint equation 



E^H Tg) = E { - ^i(^)7?72r -) - w c5(o;)^f r w + \cm^- v i1^u) 



/3 



(64) 

Since the LHS is symmetric in a/5, all terms on the RHS that are not symmetric should 
cancel. The terms from Eq. (|62j) that are proportional to A^£>^- give another constraint 



(65) 



In Eqs. (IB^jl and (}65|) . the indices a and /3 are purely perpendicular. The equation that 
defines 5"B^ is the same as Eq. (|S2>, just with the Q 13 D irac structures. 



C. Solutions to the constraint equations 

We now find solutions for the 0{\ 2 ) constraints in Eqs. ()55I59I60I64I65|) . Note that 
by careful construction of our operator basis we have ensured that each equation gives a 
constraint on a different NNLO operator. 

Eqs. (|51|55|f)4|fi5|) have implicit spinor indices, one or two vector indices, and a sum in j 
over independent structures. Since all of the equations appear between [£ n • ■ • 7i v ] they are 
only valid when the spinor indices are projected onto a 4-dimensional subspace, rather than 
the full 16-dimensional space of Dirac structures. 

It is useful to exploit the following method to determine how many independent Dirac 
structures we should have for each operator. Start by consider the three minimal structures 
that appear in the trace reduction formula, Eq. (fT3|). namely {1,75,7"}. Next for each 
case write down all possible scalar objects (v M , n M , g^, . . .) to saturate the Lorentz vector 
indices coming from derivatives in the operator and current indices, taking into account any 
symmetries. To satisfy parity and time reversal with 75, we will need to have an e-tensor, 
such as iej^j^. As long as the scalar objects are linearly independent these steps give a 
complete basis. 

At 0(X°\ a complete basis of Dirac structures for scalar, vector, and tensor heavy-to-light 
currents is |2j 

r w = i, r{! ■ =It»,v»,— }, r^ r {^ )7 ^^- 7 ¥,-nM}. (66) 

1 1 l n-v J 1 ; l. n-v n-v J 

At 0(A), there is no constraint on j( lb \ and Eq. (joTjl constrains the currents in terms 
of J(°\ To impose this constraint we need 

<^r (1) = o , ^rf 12) = o , tfrfo = g a / , (67) 
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The constraint equation causes some Dirac structures to always appear in the same combi- 
nation. We find 
1 

2 n-v 



a I a na „,q 



©Si-3) = {^7^^},|t^ + 2^ 



U (al 



( a1 - 4 ) ~ l2nV ±l ' ' ' ' 2 1X1 y± ' ' 2 



°(61-6) 



{7jr5r-4) . rfV 1 » ( 68 ) 



where T( 1 1 l/ _ 4 ) are given in Eq. which is in agreement with Ref. 17[. This basis is 

equivalent to the one in Ref. jl3]- 5 We take Q b j terms with no $ so that this choice does not 
need to be modified if we enlarge the basis for v ± ^ (see section ITllD|) . With Eq. (J58J) . 
the constraint Eq. (|51|) gives relations for the Wilson coefficients in the J^ la ' current 

^iV) = cfV), B$_ 3 (u) = C&( w ), = C&(w) . (69) 

These results agree with Refs. 0, E3 • 

At 0(X 2 ) we must solve Eqs. (|55l59l6UI64l65j) . From these equations we see that the 
currents J^ 2g \ J^ 2h \ j( 2j \ and are not constrained. The currents j( 2a \ J^ 2b \ 

j( 2c \ and j( 2d ' are all related to the leading order current J(°\ Finally the currents 
are related to the currents and J^ lb \ 

To solve the equations we will need 

5 $ CT r (1) =0, 5^=0, ^rf 2) =^, 5lT^=n° T n» (70) 

«^rg = o , a^rg = 7 ^ , ajrg = , s^rfi = + ^n"^* , 

where = — {n-v) 2 n?p. We will also need 

(71) 



°* W (61) ~~ 


a/3 f 






°* W (61,2,3) ~~ 


-#f n-v | r^ )2)3) + 7?^rf lj2)3) , 


°* W (64) 


a/3 n T 


- W (M, 2,3,4) ~~ 


a/3 ^ "pA*^ i /3 fa-p/tf 
— 9± n ' V 2 1 (1,2,3,4) + Tl * 1 (1,2,3,4) ' 


°*^(b5) 


n-v 


°*^(fe6) ~~ 


1 a/3 Lu i/l 

n-v 







where a/3 were projected onto _L directions. Note that 5*©( a .j) are easily obtained from 
these. The constraints in Eqs. (|59l6Uj) have a particularly simple solution: 

A bj (uj) = ujC'^u), A cj (u) = B aj (u), T {bj) = T {j) , T {cj) = Q {aj) . (72) 

Solutions to the other equations are slightly more involved. We present solutions to the 
constraint equations for the scalar, vector, and tensor currents in turn. 



5 Note that a structure g^n^ is redundant in 4-dimensions |29|, |31 
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1. Solutions for scalar and pseudoscalar currents at 0(X 2 ) 



The RPI constraints do not effect the allowed Dirac structures for scalar currents, so we 
have the complete sets 

T ^ = \ lT '~v) ' T(bl) = T(sl) = 1 ' T ^ = 2^ 7± ' 

_ a/3 rya/3 _ ya(3 _ _ npa/3 _ f a/3 a J3\ / 7 o\ 

L (dl) — y± J (el,2) — i (/l,2) — L (hl,2) ~ i (*l,2) ~~ l^i- >717_LJ- 

For the four quark operators, there are three possible Dirac structures in the Xn"'Xn 
bilinear, {^,^75,^7"}. In performing the matching onto SCET at a scale ~ m^, the light 
quark masses are perturbations, and for matching onto the 0(X 2 ) four quark operator we 
can set m q = 0. In this case, chirality rules out the ^7" structure which connects right and 
left handed quarks. A complete set of structures is therefore 

(T«T) aij2) = (T«T) ( ^ 2) = {^«l, ^7W}. (74) 

To solve the RPI-$ constraint, we insert the Dirac structures Eqs. (|6b17UI73|) into Eq. ()55|). 
Satisfying this constraint requires a relation on the Wilson coefficients 

4>> H = Cf >( w ) 4fH = 2uC?'{u) . (75) 

The solution for the SCET RPI-a constraint equation in (f72*|) gives 

A^{u J )=uC\ s) \u) 4?M = C\ s \uj). (76) 

To solve the SCET RPI-* constraints in Eqs. ()64I65|) . we need the additional Dirac structures 
in Eqs. (J68I71|) . On the RHS of Eq. (JSH) we observe that all structures that were not 
symmetric in af3 cancel, in agreement with the symmetry of the LHS. Solving the equations, 
the relations on the Wilson coefficients are 

J${u) = (77) 

4?K 2 ) = --cfW^ 2 ), 

002 

4?K 2 ) = -—^—Ci'\u l +u a ) - Bl?(u 1>2 ) . 
{U1+U2) 

The following Wilson coefficients of scalar currents are not determined by the RPI constraints 

4?>1, 2 ), A gl,2^1,2) , 4t 2 K2,3), Ag 2 K 2j3 ), AJJ^^Wxas). (78) 

Since the light quark in the full theory current retains its chirality in the effective theory 
current, the results for the expansion of the pseudoscalar current, q^b, are simple to extract 
from those for the scalar case, qb. The Dirac structures for pseudoscalar currents may 
be obtained by multiplying Eqs. (|73I74|) on the left by 75 and 1 <8> 7 s , respectively. The 
constraints on the Wilson coefficients of these currents are then identical. 
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2. Solutions for vector and axial-vector currents at 0(X 2 ) 



The analysis for the scalar current can be extended to the vector currents, where the 
extra Lorentz index makes ensuring that the Dirac basis is complete slightly more difficult. 
We use the method discussed in section IIII CI to count the number of terms in the Dirac 
basis prior to imposing the RPI constraints. For the case of the index a is transverse 
to v and we have 

1 : rf T n^ rf T v»} , 7 5 : {iff} , ll : R, «"} , 7^ : {n a T } , (79) 

which has seven elements. The counting for the T& cp cases are straightforward. For Y'j 
the indices af3 are _L and symmetric. We have 

1 : gfv»} , 75 : {-} , 7^ : {<??} , 7? = , (80) 

so there are four elements in the basis. Finally, for T° ^ 4 we have 

1: {gfn^tfv*}, 7 5 :{<V,i6«V}, 7 ^:{^f}, 7l : {#}, 7? : {^}, (81) 

so the basis has seven elements. 

For computations, a different basis choice is slightly more convenient. The independent 
Dirac structures appearing on the RHS of the constraint equations reduce the basis for 
Y{„i_7} by one further element. For the vector currents we find 



„ , „ ,, , rin 



T M- 6 ) = 7 M 7? , v^ T + 2g a / , —7? , {7", <A — }— , (82) 



= "I Sf {V, <>", — } , " 
1 [ n-f 2n-v 

i ; [ n-t> n-v 2n-v ) 



The index symmetrization means j^g^ = 7"$^ + 1a_9°l- I n ^q- ^82|) we have used 
Eq. (p"3|) to remove redundant structures. 

The operators j( 2a > 2b ' 2c ) bear some similarity to the complete basis of six 1/m suppressed 
heavy-to-light currents in HQET [if], 47]. The differences are due to the fact that for a 
collinear light quark we have the vector n fl available to build additional structures and from 
the fact that working in the v - 1 = frame, we do not need operators like x n iv-T>^'X r H, v . 

For the four quark operators, a basis of Dirac structures is 

( T ® T )oi-6) = ( T ® T )r,i- 6) = {— }, — 7 5 ®7 5 {7^,^,— }}• (83) 

' u- 1 o) ' (ki o) in-v n-v n-v n-v J 

Here the counting of the number of independent structures proceeds in the same way as 
for the bilinear operators, except that we start by writing down minimal structures for the 



23 



four quark operator where we impose the correct chirality on the purely collinear fermion 
bilinear. For we start with six structures, , ^75} <8> {1, 75, 7°}, and find that only the 
six terms 

^® 7 ±:{1}, ^5®7l:{<}, 1*75 ® 7s : , (84) 

are allowed, which we swap for the basis in Eq. (|83|). The analysis of discrete symmetries for 
these currents is similar to that of the four quark operators in the HQET Lagrangian |48| . 

Using Eqs. ()66I70|) . the relations for the vector current coefficients obtained by solving 
the RPI-$ constraint in Eq. (|55|) are 

4t 3 H = (&(<*) , = 2uC?'{u) , 

A^{u) = %>dP'(u>) , A®{u>) = -2CfM+2^cf V) , (85) 

The RPI-a solution in Eq. (j72|) gives 

AtUu) = coC[%(u) , 4?- 3 (u,) = C[%(u) . (86) 

Using in addition Eq. (J82j) . we find that solving Eq. (J64j) gives 

A$(u>) = -uC^\u) , A${u>) = -uC^'{u) - 2Cf (u;) , 

Ag> (u) = -uiC^Xu) + 2C#> ( W ) , < (") = (u) . (87) 

Finally, solving the second RPI-* constraint in Eq. (J65|) gives 

4i } K 2 ) = -(-)c 1 (u; 1 +lu 2 ) - B b3 (u; 1>2 ) , (88) 

\L0 2 / 

^2(^1,2) = -( — jC^wi+wa) - B m Kj) , 

4? (wi, 2 ) = -f— )c 3 (cj 1 +^ 2 ) + B b3 (u li2 ) + fi 64 (cu li2 ) , 

4? K2) = -f^— )c 1 (a; 1 +a; 2 ) + B bl (u li2 ) + 2B b3 (u lj2 ) , 

4sVi.2) = -(— ^)c 2 (wi+c<; 2 ) - 2 J B 61 (u; 1 , 2 ) - B b2 {u 1>2 ) , 
Wi+to> 2 / 

4?K 2 ) = -(^—^Csi^+uj,) - 3^3(^1,2) , 
\cg>i+u; 2 / 

4? K2) = -2B b3 (u lt2 ) + B H (u lj2 ) . 

The following Wilson coefficients of the 0{\ 2 ) vector currents are not determined by the 
RPI constraints, 

4t>i,2) > 4i-3K2) , , • (89) 

The Dirac structures for axial-vector currents which expand ^757^6 may be obtained by 
multiplying the Dirac structures in Eq. (|H2*j) by 75 on the left and in Eq. ()83|) by 1 (g>7 5 . The 
relations for their Wilson coefficients are then the same as the vector currents. 
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3. Solutions for tensor currents at 0(\ 2 ) 



The counting of the number of independent terms proceeds just as in the vector case but 
now with antisymmetric indices \w. For j( 2a \ the index a is transverse to v and there are 
ten structures 

1 : {vK^n a , gfv u] , gfn u] } , 7s : {zef n CT , tefn u] , iefVl } , 

7 [ M : {u'V, n'V, (/f 7 } , 7 1 : {vK u] } . (90) 

The bases for j( 26 ' 2 #) are simple, while for j( 2c ) we have six terms 

1 : {gl l W,gt [ W} , 75 : {ief n", ief* v*} , 7 f : {<?? } , 7 2 : {^} . (91) 

We also have six terms for 

l:{gfn^}, 7s : {ie^gf } , 7 f : {n^gf ,v^gf} , 7? : {9^, g^} , (92) 

where the identity g^e^ = —g°fe^ leaves only one term for 75 . Finally for j( 2e > 2 /> 2/l > 2 *) we 
count ten terms 

1 : {gfnK"\g^gf} , 7s : {ie^gf, iefnK^} , 7 f : {n^gf, v^gf} , 
ll ■ {9?n»\g^} , ^ : {g^n v \ gf^} . (93) 

Again only j( 2a ) has its basis of Dirac structures further restricted by the RPI-$ constraint 
in Eq. (I55j) . which reduces the basis by two terms (since only eight linearly independent 
Wilson coefficients appear in Eq. ()96|) below). For the complete set of Dirac structures for 
tensor currents we find 



1 (ol- 



n-v n-v 



-p/ii/ _T_\ ly^tiv fpM^ 

^-^n-vS ' (6>9l— *) ~ (1-4) 



1 (el-10) 



(^±1(1-4) > 7±7±1 (1 _ 4 ) , 2n .^Ji 7 , 2n-v 7±9± J ' 

T Q/3 ^ - lla^T^ ryO-rflY^ V 3 o a[A V] V 3 o a[/ V 1 l C941 
1 (/,/i,a-lo) ~~ (1-4) ' 7_l7± j - (i_4) j 7±i/± 7 ) 7j_i/j_ u J) l y/± J 

where = y a/3 — u a ?A Similarly, for the tensor four quark operator currents, a complete 
basis is 

(T0T)-_ 1O) = (T«T)-_ 1O) 

J- ® rr 4) , ^-T 5 8 7 5 rS1 4 , , io'" ® 1 , i<r<"" 75 ® 75) , (95) 
n-v K ' n-v K ' J 



25 



where just as for the vector case we have made use of chirality. 

The relations for tensor Wilson coefficients obtained by solving the RPI-$ constraint 
equation are: 

4?( W ) = C 1 H, 4?H = C 2 {u) , (96) 

A®(u>) = C s {u>), A®(u>) = C A (u>), 

A®{u) = 2uC[(u) , A<2(cu) = 2uC' 2 (u) , 

4?M = -2C 3 (oo) + 2uC' z {oo) , A®{u) = -2C 4 (uj) + 2ujC' 4 (uj) . 

The relations for Wilson coefficients from the RPI-a constraint equations are 

41 4 M = wCS(a;), 41» = C^M- (97) 

Finally, solving the RPI-* constraint in Eq. (JH%|) for the tensor case gives 

4? (co) = -uC?'(u) , Ag {u) = -uC®\u) - 2Cf ( W ) , 

4g(u;) = -uCf^) + 2Cf (w) , A®(u) = -uCf'(u) + 2Cf 



(«) = -2Cf ( W ) , («) = 2Cf ( W ) . (98) 

while the constraint in Eq. (|65jl has the solution 

4? K2) = -(-)c?W<"a) - 2S«(^) , 
!(*)/-,. \_ f m \Mt) f \ D {t), 



42 H2) = -l — )C™(u 1 +u 2 )-B™(u 1 ri+B K (u li2 ) 



2) , 



4? ("1,2) = -(-)cl\u 1+ u 2 ) - B®( Ulfl ) - K 

42 K2) = -(-)C^ (o^+wa) - 2S«(^) + 4W) + B mM > 
42 H*) = -(^McfW^ 2 ) - Bg>( Wl>2 ) + 4i$K 2 ) , 
48 Ha) = - (-^-)cf ( Wl +a*) - 2if (u; 1)2 ) + Bg^) + 2<> (0,1,2) , 
4? H2) = -(-^)cf (wi+a*) + 3^?K 2 ) , 

42 H2) = "(-?-) C4' } ("i+wa) + ^SW " 32$ HO , 
A2(a; 1>2 ) = -2Sg(a; lj2 ) + ^f(a; 1)2 ), 

4?oH0 = 45^(^,2) - 25$ (0,1,2) + B^M . (99) 

The following Wilson coefficients of the 0(X 2 ) tensor currents are not determined by the 
RPI constraints 

4I10K2) , 414(^1,2,3) , 4,Lio(^, 2 ) , 4LioK2, 3 ) . (100) 



26 



D. Absence of supplementary projected operators at C(A 2 ) 



Here we show that the analysis above on the surface v ± = is complete by showing that 
there are no supplementary projected operators as defined in section lll CI The analysis of the 
proceeding section makes this simpler, since a complete set of relations have been derived for 
all currents Jj 2b ~ 2e \ Thus, we only need to worry about supplementary projected operators 

generated by transforming the currents jj 2a ' 2 ^ 2k \ To simplify our proof we first swap all 
factors of n-v for l/(n-v). 

First consider the SCET RPI-II transformation at O(X ) for these currents. At this 
order we have 

n^^n^ + e^, 7^ -»• 7^ - — f ± - — i , -> — - e ± -v ± - — n-v , 

n-^-y^i. S^^-yei-Bx. (101) 

We use the convention where all indices a (3 are _L for the field structures and Dirac structures 
in Jj 2 ^~ 2k \ Now due to the contractions of the a and (3 indices only the transformations 
on and 7^ can contribute for these operators (there are no n's or _L's in the j( 2a ) case). 
The transformation related to their labels Ui is 0(A) and need not be considered and the 
field transformations cancel. Thus, the only terms that appear in an RPI-II relation are 
those whose Dirac structure transforms, 8i{Ti a j, g ,h,i,j,k) 0- However, with our choice of the 
complete basis of Dirac structures on the v± = surface, the structures for these currents 
all have zero transformations. In this regard it was important to take a basis with no factors 
of f/i. Away from this surface we must add to our basis of Dirac structure by including 
additional v± dependent terms and it is only these terms that can have additional relations. 
For example, factors of v± are induced when we reduce a basis that includes factors of ^ 
using the trace formula in Eq. (|13p. The same is true with our choice of the basis of 
currents. 

Finally consider whether the transformations RPI-* and RPI-$ induce SPO's or equiva- 
lently SCET RPI-I and HQET RPI. Since the n transformation in RPI-$ did not enter at the 
order we are working it is apparent that there are no SPO's from the HQET RPI. Examining 
the results of the RPI-I transformations we find that none of the jj 2a ' 2 ^ _2fe ) currents have 
O(X ) transformations (since the Dirac structures transform at O(X) and the field structures 
that do transform all cancel out). 

Thus the results derived in the previous section give the complete set of RPI relations 
for the 0(X 2 ) currents when v± = 0. 



IV. CHANGE OF BASIS AND COMPARISON WITH TREE LEVEL RESULTS 



In expanding the heavy-to-light currents, two different bases of operators are useful. At 
tree level it is convenient to write the result for the currents in terms of collinear covariant 
derivatives, giving one basis. For the derivation of RPI relations and factorization theorems, 
a basis such as the one in Eq. (|%nj) is more useful. 

The tree level matching of the full theory current qTh onto SCET currents was done to 
subsubleading order in Ref. In deriving Feynman rules we find the momentum space 
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version more convenient so we use the equivalent result from 0] 



-WO) , --(la) , --(16) , --(2a) , --(26) , --(2c) , --(2d) , --(2d) , --(2/) 



^rV 6 -> j (U > +T a> + T"> +T za> +j 



x 



x 



x 



-X 



(102) 



where 



T (0), 
J [UJ 

J (a; 



-f(2a) , 
J [UJ 



-y(2c), 
J (a; 

T (2e) , 

J (^1,2 

-7(2/), 



1 



(XniX>^ A _) uj Q(a)a'Hv 



^ r Y > '-T-) 1 a 'Ti 



-f(16)/ \ 

J (^1,2) 

F(2&), 



-1 

m 
1 



(103) 



UJ 



n-v ign-B c \ — 



V n-v 



J (^1,2) 



ign-B c \ — 



-1 _ 

m \ n-v / uj 2 



— fc^cl) wl (^fl).J(e)^ ^ , 



f(0) . 



The r in J is simply the Dirac structure of the full theory current. The Dirac structures 
that appear in the subleading currents are 



©(a), 



lot o- ' 

In ■ v 



6 



(6) a 



2n ■ v 



To 



(104) 



and 



T( a )a = T7, 

t w = r 



T 
a ) 



4 



I (6) a = 7a 7T 1 ' 

In ■ v 

T( e ) a/3 = 7a g 1 2 7 ^ ' 



T 



(<0 



T(/) a/3 = T7 Q L 7^ 



(105) 



Each of the operators J has unit Wilson coefficient at tree level. By re-expressing these 
operators in the basis of operators presented in this paper, we determine the tree-level 
Wilson coefficients of our currents. This provides a check of the RPI relations. 



A. Conversion 

In terms of our basis, the leading order tree level current is given by 

J iuj J (0) H = Jdoo J{ 0) (u) . (106) 

This result holds for all five Lorentz types, T = {1 , 7 s , 7^ , 7 5 7 M , ia^}. For the remainder of 
this section, we will suppress the explicit cu-dependence of our basis J(o»)'s as well as the ap- 
propriate integrals /[do;*] whenever results hold equally well as integrals or as densities. For 

example, Eq. ()106|) would be written simply as = If the Lorentz type (s,p,v,a,t) 
of the current is not specified, the same result holds for all five types as above. 
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For the 0{\) currents, the relations differ for the scalar, vector, and tensor cases, 

m -jMiuuUt), (107) 

[j} lb V^2)-2jfVi^2)], 

J J uj\-tuj 2 

Jdui jf a \uji) = jduji J^ a \ui) + fdu 1>2 1 1 ™ t t [-Ji lb \u 1 ,u 2 ) - 2J 2 (lb) (o;i, u 2 )] , 



JdUi J^M = Jdui Ji a) {u)i) - Jdu 1:2 
dui (ui) + Jdco 1;2 
dui J^ a \ui) + du lj2 



0Jl+UJ 2 
U}i+U 2 

rn 



J 



(ib) 

s,p 



0, 



J 



(16) 
v,a 



-J, 



(16) 



3 > 



LO\+UJ 2 

j! 16) 



2 + 4 ib) + 2j[ ib) + 2 4 ib) 



The last line of relations are true as integrals or as densities. For example Jv]a\oJi j2 ) = 
— J^ lb \uji )2 ). At 0{\ 2 ) the relations between the two forms of subleading currents are the 



-j(2a,2b,2c) 

same tor all J currents 



j(2o) = j(2a) 



7 (26) = jf\ fal™ = J< 



duii 



m 



UJ 2 



_jf ) + , (108) 



f( 2c ), 



where in the last relation the arguments of J^ J (ui) and j[ 2e,29 \uji j2 ) are implicit. The 
remaining currents come in different combinations depending on the Dirac structure. For 



J we have 



s,p 



J 



(2d) 



j(2e)_ j(2g) 
J 3 J 3 i 



— (2e) 

For the scalar J currents, 
J du h2 J ( *p =J dw K 



j(2e)_ j(2e) 



+ duj h2 



J 



f 



(2d) 



-4 2e) +4 29) 



UJ 2 + UJ 3 
LUX+UJ2 



h 



(109) 



, (110) 



with similar relations for the vector and tensor cases (suppressing the integrals for conve- 
nience) , 



7 (2e) _ 
J v,a 



J 



uj 2 +u; 3 
uj 1 +u 2 



(2c) T (2e) . T (2e) 



(2e) 



+ 



(2f) _ j(2f) _ T (2f) , T (2f) 



j(2h) j(2h) j(2h) j(2h) „ j(2h) 



J 



(2e) 



jj 2e) - 4 2e) + jf e) + 4 2e) 



+ 



./ 



+ 



UJ 2 + U03 



jr-"+Jk 



h 



(2/) 



Finally for J (2/) , 



(2/) = 

s,p 

(2/) 
v,a 

j! 2/) = 









e — > /i 


0J 2 









111) 



(112) 



-J. 



(2e) 7 (2e) 



(2e)l _ ^1 
6 -I 0J 2 

- Jf e) - 4 2e) + jf e) + 2jf e) 



oo 2 



e — > /i 



/i 



where the suppressed integrals are the same as for J 



T(2e) 
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B. Wilson coefficients at tree level 



Inserting Eqs. (jl06H112j) into Eq. (jl02|) . we can read off the tree level Wilson coefficients 
of our basis. For example, since J is the only term at leading order we have C[ d \uj) = 1 
and Cj%(u) = for d = s,p, v, a, t. 

For scalar currents, the non-vanishing tree-level Wilson coefficients are 

C[ s \ui) = 1 , B$(u) = 1 , Bg(uv) = , (113) 

and at 0(X 2 ) 

4?k 2 ) = --, 4iW = i + -, 

UJ 2 1 UJ 2 

1 <JJ 2 y L) 2 

As), v ^l + 2cJ2 + ^3 As), v + + 

^{(^1,2) = , A\ 2 \u la ,3) = • (114) 

LOi+L0 2 UJ\-\-UJ 2 

The same results hold for the pseudoscalar currents. To C(A 2 ), the values of the Wilson 
coefficients for vector currents that do not vanish at tree-level are 

b^M = , sffKO = -=^- , s&W) = -1 , (H5) 



and 



u 2 



A%(u 1>2 ) = -2, < j K 2 ) = 3, 

4?K2)=2, ^(a;^) = 1, 

4?(a; 1)2 ) = -1 + ^ , A$(u lj2 ) = -1 + ^ 

4; } k 2 ) = i-— , 4?k 2 ) = -2, 
^2 

,(«)/ ^2+^3 ^1-^3 



( „) a;i-a;3 (v) _ -2^-^+^ 

Am (Wl ' 2 ' 3) " ' Afc6 (wi - 2 ' 3) " Wl +u; 2 ' 

A«( Ww ) = -2^. (116) 
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The same results hold for the axial vector currents. Finally, for the 0(X 2 ) tensor currents 
we have nonvanishing coefficients 



cfV) = i, B£M 



(*)/ 



2 - 



m 



-2m 



B® M = l, B W (o; 1)2 ) = 2 , B W = 2 



i(t) 



(*)/ 



117) 



and 



j(t) 

,(*) 



.4 
.4 



At) 
At) 

A gl 

At) 
A h\ 

(*) 
h5 



A 



u) = 1 
^1,2) = 

vi?) = 
^1,2) = 

^1,2) = 
^1,2) = 
^1,2,3) 
^1,2,3) 



-f 

f - 

m 



m 

U> 2 



UJ 2 



<jj<l 

UJ1+UJ2 ' 
— Cl>!+Cl> 3 

UJ1+UJ2 



4?M = 

4?K 2 ) 
4*3(^-2) 

4t(^, 2 ) 
41(^i 
4?K 



2wi 

Co>2 



2,3 J 



2,3J 



^1-^3 

Wi+u;2 ' 

— 2^1—^2 +CO3 

UJ1+UJ2 



(118) 



It is straightforward to check that these results all satisfy the RPI relations from section lTlI CI 
providing a cross-check on those results. 

C. One-Loop Results 



The relations from section 1111 CI apply at any order in perturbation theory, so they can 
also be used to determine one-loop values_for certain coefficients. For the LO currents the 
one-loop coefficients in MS at fi — m are 



Cl v \u 



a s (m)Cp [ „, „ T . . AV 21n(a)) 7r 2 

1- sV , ^ <!21n 2 (t2>) +2Li 2 (l-cu) - ; v / + — , . 



47T 

a s (m)CF 
a s (m)CF 



2ln 2 (^) + 2Li 2 (l-tD) + ln(w) 
2ln 2 (t2>) + 2Li 2 (l -w) + ln(w) 



1-w 12 

3^ — 2 



47T 

a s (m)CV 

4^ ' (1 — cD) 2 J ' 

a s (m)Cjr f (1 — 2a>)u) In (a)) a) 



1 -a; 
4^-2 



7T 

12 

^ 2 



6 , 



l-^ + 12 +6 >' 



2 ^ 2cD ln(cD) 



4tt [ (l-^) 2 
a s (m)Cjr f— 2a>ln(a>) 



47T 



1 — cD 



1 - a- 1 ' 

cf(a) = o, 



(119) 
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where cj = cj/m and Cp = 4/3 for color SU(3). The quark-gluon-antiquark operators 
J^ 16 ' have coefficients that are not fixed by RPI, and these were determined by a one-loop 
matching in [29, 30]. Thus all 0(X 0,1 ) currents are known at one-loop order. The expressions 
are fairly lengthy, and so we do not repeat them here. Using their results and our Eqs. ([751 
1771) . JEMEEJ), and flUED , the coefficients of the currents j^^cM,^) are algo determined 
at one-loop order. 

We give the scalar current case as an example. For the scalar current, the coefficient at 



m is 



£(bl)(&l,2, 1) 



-^^[-(ln^-ln^-lnf^ 

CJ 4:TT \_UJ2 { 

In cj In o>2 
1— cj 1 — cj 2 



2 



0>i/ 

u) 2 hiu) 2 2(1— Ui) 

(1-CJ 2 ) 2 Cj\Cj2 



CJ 



1-UJ 

{Li 2 (l-uV)- 



lncj 
Li 2 (l-o>i)} 



cjicj 2 

47T 
1 — 0>i 



r . 7T 2 4 1 

|Li 2 (l-^)--)- s - — 



CJ 2 



-1 



— <^ ln^-ln^i-ln — +— ln — + 



CJlCJ 2 



CJ 2 

{Li 2 (l- 



( u 



In d> 2 

1-CJ 2 
2 . 



-W) 



Li 2 (l-wi)} - ^^{Li 2 (l-^ 2 ) - 



47T CJ 



21n 2 (cj) + 2Li 2 (l-cj) 



cjicj 2 

1 



7T 

^ f 12 
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where Ca = 3, cji )2 = cji >2 /m, = cDi + cj 2 , and we have transformed to our basis. We will 
also need the derivative of the LO scalar currents coefficient 



gscj 



a s (m)Ci 
27rm 



•1 + cj + (2 - 4cj + cj 2 ) lntj 
Qft -Cj) 2 



121) 



Now in section IIIC we derived the following results for the 0(X 2 ) currents 

4?M = cfM, Agfa) = %,(#»'> 



4?H 



ujCI 



(«)', 



CJ J 



x a2 



CJ 



(122) 



^V) = cf(cj), = -cjCT'M 



(*)/ 



.(«)'( 



n? 



c; s) (cj 1+ cj 2 ), ^(cj 1i2 ) 



lJ 2 



??? 



C; S) (CJ!+CJ 2 )- B^Wy). 



(cji+cj 2 ) 



Combined with Eqs. ()l 191112 1|) . these relations determine the coefficients at one-loop order. 
The results for the jf a 2e ^ vector and tensor currents at one-loop order are easily obtained 



in the same manner. 



V. CONCLUSION 



In this paper we derived a complete basis of scalar, vector, and tensor heavy-to-light 
currents at next-to-next-to-leading order in the power counting, 0(X 2 ). Building on the 
approach in Ref. [il| where one takes v± = from the start, we constructed the full set 
of RPI relations that leave us on this surface. The completeness of deriving RPI relations 
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projected on a surface was analyzed. With a careful choice of Dirac structures in our analysis 
of heavy-to-light currents at 0(X 2 ) it was demonstrated that the projected RPI gives the 
full set of constraints. We also investigated the path dependence of Wilson lines in order to 
clarify what conditions they must obey to give the correct cancellation of usoft gluon effects, 
and to demonstrate the manner in which results are independent of the choice of boundary 
condition. 

A simple method for counting the number of Dirac structures in the basis for any operator 
with d = 4 was given. 6 Several types of reparameterization invariance provide restrictions on 
the structure of these currents, which we formulated as constraint equations on the allowed 
Dirac structures and Wilson coefficients as given in Eqs. (f5Tj) . (JH5|) . (jSHJ), (ffiO|. and 
(Ifj5jl . We expect that a similar setup with constraint equations and projected surfaces will 
be useful in deriving RPI relations at higher orders in A and in deriving results for non 
heavy-to-light currents. 

Our main results are contained in the solution of the constraint equations as given in 
Eqs. (J7H77D, ijBSISBjl . and (jHHEHD- These results determine the coefficients of five of the 
eleven NNLO operators, j( 2a ' 2b < 2c > 2d > 2e ^ f or various Dirac structures indicated by j and at 
any order in perturbation theory, in terms of the coefficients of NLO and LO operators. This 
determines 7, 23, and 32 Wilson coefficients for the scalar, vector, and tensor heavy-to-light 
currents respectively. Results at tree- level and one- loop order were discussed in sections TlVBI 
and llVH Finally, the operators j( 2 f?9?h,2i,23,2k) defined in Eqs . (j4Q|4ij) together with the 
Dirac structures in Eqs. (|73l74l82l83l94l95j) were shown to not be constrained by reparam- 
eterization invariance. 
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